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ABSTRACT 

We present multi-frequency observations of a sample of 15 radio-emitting Broad Ab- 
sorption Line Quasars (BAL QSOs), covering a spectral range between 74 MHz and 
43 GHz. They display mostly convex radio spectra which typically peak at about 
1-5 GHz (in the observer's rest- frame), flatten at MHz frequencies, probably due to 
synchrotron self-absorption, and become steeper at high frequencies, i.e., v >20 GHz. 
VLA 22-GHz maps (HPBW ~ 80 mas) show unresolved or very compact sources, with 
linear projected sizes of ^1 kpc. About 2/3 of the sample look unpolarised or weakly 
polarised at 8.4 GHz, frequency in which reasonable upper limits could be obtained 
for polarised intensity. Statistical comparisons have been made between the spectral 
index distributions of samples of BAL and non-BAL QSOs, both in the observed and 
the rest-frame, finding steeper spectra among non-BAL QSOs. However constraining 
this comparison to compact sources results in no significant differences between both 
distributions. This comparison is consistent with BAL QSOs not being oriented along 
a particular line of sight. In addition, our analysis of the spectral shape, variability 
and polarisation properties shows that radio BAL QSOs share several properties com- 
mon to young radio sources like Compact Steep Spectrum (CSS) or Gigahertz-Peaked 
Spectrum (GPS) sources. 

Key words: quasars: absorption lines — radio continuum: galaxies — Polarisation 



1 INTRODUCTION 

About 15 per cent of optically-selected quasars exhibit 
Broad Absorption Lines (BALs) in the blue wings of the 
UV re sonance lines, due to g as with outflow velocities up to 
0.2 c (|Hewett fc Foltdl2003h . A classification of BALs has 
been done according to their association with highly ionised 
species like Civ, Nv, Siiv (HiBALs) or, in addition, with 
low-ionisation species like Mgn or Aim (LoBALs). These 
broad troughs are the most evident manifestations of quasar 
outflows but many other intrinsic absorbers such as Associ- 
ated Absorption Lines (AALs), mini-BALs or high velocity 
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Narrow Absorption Lines (NALs) are also related to them 
l|Hamann fc Sabrall20oJ : ICangulv fc Brotherton]|2008l ) . 

The study of BAL quasars (BAL QSOs) has become 
increasingly important in the last years. Outflows are not 
only useful as probes of the physics in the AGN environment, 
but also because they might be an important piece in many 
other puzzles like enrichment of the intergalactic medium, 
galaxy formation, evolution of the AGN and the host galaxy, 
cluster cooling flows, magnetisation of cluster and galactic 
gas and the luminosity function of quasars. 

At the moment it is not clear how many quasars host 
outflows and why these can only be seen in a fraction of the 
quasar population. The most popular hypotheses proposed 
to explain this mainly differ in the role given to orientation. 
The main premise in the 'orientation scenario' is that the 
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BAL phenomenon might be present in all quasars but in- 
tercept ed by only ^15 per ce nt of the lines of sight to the 
quasar l)Wevmann et al.lll99lf ) , e.g. within the walls of a bi- 
funnel centred on the nucleus (|Elvdl200Ch . This scenario is 
supported because in most respects, apart from the BALs 
themselves, BAL QSOs look similar to normal quasars. The 
small diff erences from non-BAL QSOs, e.g. slightly redder 
continua peichard et al.ll2003bl ) and higher optical polari- 
sation, could be a consequence of a preferred orientation of 
viewing angle. 

The alternative hypothesis, the so-called 'unification by 
time', assumes that BALs appear during one or maybe more 
short periods in the lifetime of quasars. The observed rate of 
BAL QSOs would then be a function of the percentage of the 
quasar lifetime in which these outflows show up. This evolu- 
tionary scenario has been discussed bv lBecker et al.l (|2000l ) 
on the basis of the radio properties of the ir Faint Images Ra - 
dio Sky at Twenty-centimeters (FIRST. [Becker et al.lll995l ) 
BAL QSO sample. They note that 80 per cent of their BAL 
QSOs are unresolved at a 0.2-arcsec scale, with both flat 
and steep spectra. This latter group resembles the Compact 
Steep-Spectrum Sources (CSS) and their Gigahertz Peaked- 
Spectrum (GPS) sub-class. As a consequence, they suggest 
the picture in which BAL QSOs represent an early stage in 
the development of quasars. 

This last scenario finds motivation on the collection 
of theoretical and phenomenological works that in recent 
years have contributed to build a first consistent picture 
of early AGN evolut ion. As an example, recent simulations 
popkins et al.ll2005l ) describe how quasar evolution can be 
understood as the result of merger events, where the super- 
massive black hole grows via accretion and, when it finally 
becomes an active 'protoquasar', expels out the surrounding 
material deposited by the merger through powerful winds. 
The fingerprints of these winds have been manifested not 
only as BAL troughs in UV spectra, but also as emission line 
outflo ws and sometim es associated to HI absorptions (see 
e.g.. iHolt et al. 2008). To that respect, even if the fraction 
of optically-selected quasars hosting BAL systems is esti- 
mated to be about 10-20 per cent of the quasar population, 
it has been argued that these samples might be substan- 
tially biased against BAL QSOs, as suggested by larger frac- 
tions of BAL Q SOs found a mong infrared-selected quasars 
|Dai et al.ll20of ). In addition. [Ganguly fc Brothertonl (|2008l ) 
have recently estimated that ~60 per cent of AGNs could 
host outflows when the different kinds of intrinsic absorp- 
tions are considered. 

Radio emission has become, in fact, an important 
additional diagnostic tool when studying the orientation 
and evolutionary status of BAL QSOs. For a long time 
radio-loud BAL QSOs were believed to be e xtremely rare 
amon g the population of luminous quasars (|Stocke et al.l 
1 19921 ). With the advent of large comprehensive radio sur- 
veys it has however become clear that radio-loud BAL 
QSOs are common. However the fraction of BAL QSOs 
seem to vary inversely with the radio-loudness parameter, 
with the radio-bri ghtest quasar s bein g 4 times les s likel y 
to exhibit BALs Pecker et all l200ll ; iGregg et all l200fj) . 
The evolutionary picture is, in fact, supported by some re- 
cent works coming from radio observations. For instance, 
EVN observa tions at 1.6 GHz of a few BAL QSOs from 
Becker's list l| Jiang fc Wandl2003[ ) show that they present 



a still compact structure at these resolutions, with sub- 
kpc projected sizes, and a variety of orientations accord- 
ing to their radio geometry. In addition, the radio power- 
ful CSS BAL QSO 1045+352 has also been observed with 
MERLIN and VLB A, and its complex radio structure has 
been interpreted as evidence of radio-i ntermittent activity 
l|Kunert-Bairaszewska fc^M arecki 2 0071 ). 

Some other works address the question o f the orien- 
tatio n of BAL QSOs through radio variability (IZhou et al.l 
120061 ; iGhosh fc Punslvl I2007T ) . They compare the 1.4-GHz 
flux densities of all quasars in t he Sloan Digital Sky - 
Survey (SPSS) quasar catalogues ([Schneider et all 120051 ; 
ISchneider et all 120071 ) with detections in both the FIRST 
and NRAO VLA Sky Survey (NVSS. ICondon et al.|[l998h 
surveys. The most strongly variable quasars have probably 
their jets closely aligned with the line of sight, being the 
relativistic beaming responsible for such high flux density 
variations. Some of those variable quasars present in addi- 
tion the BAL phenomenon which is considered as a proof 
for the existence of at least a sub-population of polar BAL 
QSOs. At the same time, a few examples have been found of 
BAL QSOs asso ciated to radio sou rces with an extended FR 
II morphology ([Gregg et alj|2006h which suggest for those 
the opposite edge-on orientation. 

This paper reports the first results of a study in which a 
systematic approach is taken to characterise the radio-loud 
BAL QSO population. The final goal is to locate these ob- 
jects in the framework of an adequate orientation and/or 
evolutionary status scenario. In order to do this, the ra- 
dio spectra and polarisation properties of a small sample of 
radio-loud BAL QSOs will be presented. 

In Section [2] some samples of radio emitting QSOs ex- 
isting in the literature are summarised, and among them a 
small sample of BAL QSOs is selected for multi-frequency 
continuum and polarisation observations. In Section [3] the 
observations, data reduction and the treatment of the er- 
rors are described. In Section [3] the results of the observa- 
tions are presented. The radio morphologies of BAL QSOs, 
their polarisation properties and the variability of some of 
these objects with 2-epoch observations are analysed. In ad- 
dition, the shape of the radio spectra is described and the 
radio spectral index distributions of radio quasars with and 
without associated BALs are compared. All these results 
are reviewed together in Section [5] and compared with other 
findings trying to put emphasis on those aspects related to 
the orientation vs. evolutionary dilemma. Some of the lim- 
itations of the analysis are noted and solutions to some of 
these caveats are proposed. Finally Section [6] summarises 
the main conclusions. 

The cosmology adopted within the paper assumes a 
flat universe and the following parameters: Hq=7Q km s _1 
Mpc -1 , S7a=0.7, Qm~0.3. The adopted convention for the 
spectral index definition is S v oc v a , except where otherwise 
specified. 



2 SAMPLES OF BAL AND NON-BAL QSOS 
USED 

The first s ignificant sample of radio BAL QSOs was pre- 
sented by iBecker et all (|200fj| ). which included 29 BAL 
QSOs. These were optical identifications of FIRST radio 
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Table 1. Sample of 15 radio- loud BAL QSOs studied in this paper. 



ID 


RA 
(J2000) 


DEC 
(J2000) 


r ro 
(arcsec) 


z 


^peafc 

(mjy/beam) 


E 


log(L 5 GHz) 
(W Hz- 1 ) 


Ref 


Type 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


0039-00 


00:39:23.18 


-00:14:52.6 


0.22 


2.233 


21.2 


19.50 


26.33 


1 


HiBAL 


0135-02 


01:35:15.22 


-02:13:49.0 


0.31 


1.820 


22.4 


16.79 


26.22 


3 


LoBAL 


0256-01 


02:56:25.56 


-01:19:12.1 


0.66 


2.491 


27.6 


18.57 


26.54 


3 


HiBAL 


0728+40 


07:28:31.64 


+40:26:16.0 


0.33 


0.656 


17.0 


15.27 


24.97 


2 


LoBAL 


0837+36 


08:37:49.59 


+36:41:45.4 


0.20 


3.416 


25.5 


19.16 


26.66 


6 


LoBAL 


0957+23 


09:57:07.37 


+23:56:25.2 


0.12 


1.995 


136.1 


17.64 


27.05 


2 


HiBAL 


1053-00 


10:53:52.86 


-00:58:52.7 


0.14 


1.550 


24.7 


18.03 


26.02 


4 


LoBAL 


1159+01 


11:59:44.82 


+01:12:06.9 


0.07 


1.989 


268.4 


17.30 


27.31 


1 


HiBAL 


1213+01 


12:13:23.94 


+01:04:14.7 


0.13 


2.836 


21.5 


19.69 


26.65 


4 


HiBAL 


1228-01 


12:28:48.21 


-01:04:14.5 


0.27 


2.653 


29.4 


17.74 


26.62 


1 


HiBAL 


1312+23 


13:12:13.57 


+23:19:58.6 


0.10 


1.508 


43.3 


17.13 


26.32 


2 


HiBAL 


1413+42 


14:13:34.38 


+42:12:01.7 


0.10 


2.810 


17.8 


17.63 


26.42 


2 


HiBAL 


1603+30 


16:03:54.15 


+30:02:08.6 


0.28 


2.028 


53.7 


17.61 


26.61 


2 


HiBAL 


1624+37 


16:24:53.47 


+37:58:06.6 


0.05 


3.377 


56.1 


17.62 


27.06 


5 


HiBAL 


1625+48 


16:25:59.90 


+48:58:17.5 


0.02 


2.724 


25.3 


17.41 


26.59 


6 


HiBAL 



Columns are: (1) source ID; (2,3) Optical coordinates in J2000.0; (4) angular distance between the optical and FIRST 
radio positions; (5) optical spectroscopic redshift; (6) FIRST 1.4-GHz flux density; (7) APS E magnitude corrected 
from galactic extinction; (8) 5-GHz rest-frame luminosity, computed using the radio spectra presented in this paper; 
(9) reference wher e the BAL QSO app ear s in the literature an d ( 10) BAL class i ficatio n a s given in that referenc e. 
References key: 1 - [Menou et al.1 [12001!) ; 2- lBecker et alj 1120001) ; 3- IBecker et alj fcOOiT) ; 4- lReichard et alj ll2003al) ; 5- 
iHolt et al.l l|2004h : 6- Identified directly from SDSS-DR3 database. 

sources in the FIRST Bright Quasar Survey ( FBQS) show- 
ing B AL features in their optical spectra. iBecker et al.l 
(200CJ) performed radio observations of this list at two fre- 
quencies, 3.6 and 21 cm, in order to determine the radio 
spectral index of these sources. 

One of the aims of our work is to present a represen- 
tative sample of radio-loud BAL QSOs with radio spectra 
covering a wide range in frequency (i.e. between 74 MHz and 
43 GHz). The selection was based on the available samples of 
radio-loud BAL QSOs in the literature when the project was 
started. The main selection criterion was a cut in flux den- 
sity of Sia GHz > 15 mjy in order to facilitate total power 
detections at high frequencies and also the detection of po- 
larised emission from those sources with a higher degree of 
linear polarisation. 

We have defined a sample of radio-loud BAL QSOs 
[RBQ sample] comp osed of 15 sources. Seven of them were 
present in the list of IBecker et all d2000h or its extension in 
the southern galactic cap (| Becker et al.ll2001n . The remain- 
ing eight objects are SDSS quasars classified as BAL QSOs 
by visual inspection and associated to FIRST s ources within 
a 2" r adius. Two of t hese were presented by iMenou et al.l 
|200ll ). three more by iReichard et all (|2003al ) and the un- 
usual BAL QSO 1624+3 7 in the sample was discovered in a 
survey for z~4 quasars (|Holt et al.ll2004l ). These 13 objects 
are the only BAL QSOs present in the mentioned samples 
with 5*1.4 ghz > 15 mjy. Finally, the other two BAL QSOs 
satisfying this flux density criterion were identified by one of 
us directly from the third data release of the SDSS (DR3). 
Table Q] shows the main properties of the RBQ sample. Fig- 
ure [T] shows flux density at 1.4 GHz versus redshift for all 
BAL QSOs belonging to the mentioned samples. Sources 
above the horizontal line of 15 mjy define the RBQ sample. 

At a later stage, iTrump et al.1 l|2006l ) presented a cata- 
logue with 4784 SDSS BAL QSOs selected using a more 
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Becker et at (2000) 
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Figure 1. Distribution of known radio-loud BAL QSOs in log S 
and z before going on-sky. The sample to be studied in this paper 
was selected by taking all those with Si ,4GHz > 15 m Jy- 



homogeneous selection criteria. They used an automated 
algorithm to fit the continuum and measure the Al (Ab- 
sorption Index, lHall et all |2002|) in all quasars from the 
SDSS-DR3 quasar catalogue (|Schneider et aill2005l ). In fact, 
all SDSS quasars in the RBQ sample belong to the list of 
ITrump et ail l|2006l ). except 0039-00 which is not included 
in the SDSS-DR3 quasar catalogue, but appears in the mor e 
recent SDSS-DR5 quasar catalogue l|Schneider et al.ll2007l ). 
However, some care has to be taken when using the list of 
ITrump et al.l ll2006l). Accord i ng to the more subjective clas- 
sification by iGangulv et all ((2007), the absorption present 
in about 15 per cent of the quasars classified as BAL QSOs 
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Table 2. Summary of the observations. 



4.0 
3.5 
s 3.0 

1-3 

e 



X Becker etal. (2000) 

O RBQ sample 

* Vigotti et al. (1997) 



- • * # #- * 

. ***•». * 




Run 


Date 


Telescope 


Frequencies (GHz) 


1 


26-30 Jan 05 


Effelsberg 


4.85, 8.35, 10.5 


2 


15-16 Jul 05 


Effelsberg 


4.85, 8.35, 10.5 


3 


27 Oct 05 


Effelsberg 


4.85, 8.35, 10.5 


4 


23-26 Jun 06 


Effelsberg 


2.65 


5 


13 Feb 06 


VLA(A) 


8.45, 14.5, 22.5, 43.5 


6 


7 Mar 06 


VLA(A) 


8.45, 14.5, 22.5, 43.5 



the redshift range between 0.096 and 2.757. Figure [5] also 
illustrates how the RBQ sample represent intermediate flux 
densities, between the flu x density coverage s panned by the 
relatively faint sample of (|Becker et al]|2000|) and that cov - 
ered by the list of brighter quasars from lVigotti et all l|l997h . 



Figure 2. Same plot as Figure [T] in which the sample of normal 
radio-loud quasars from Vigotti et al. (1997) is shown, together 
with the samples of radio-fainter BAL QSOs. 



in Trump's list might actually be due to intervening sys- 
tems instead of having an intrinsic origin. So me examples of 
these false positives are shown in Figure 2 of iGangulv et all 
l|2007f ). 

In this work the radio spectra of BAL QSOs will be 
described, and in particular, the spectral index distribution 
will be compared with that of a non-BAL QSO sample. This 
is an interesting exercise because the spectral index is known 
to be an statistital indicator of the orientation of quasars, 
and this comparison can give some information about the 
orientation of BAL QSOs. This comparison sample of non- 
BAL QSOs has b een extracted from the complete B3-VLA 
quasar catalog ue jVigotti et al.lll997h . The B3-VLA survey 
l| Vigotti et al.| [T989) consists of 1049 sources selected from 



the B3 survey (|Ficarra et allll985l ) in five flux density bins 
between declination 37 and 47 deg. It was designed to be 
complete down to flux densities of 100 mjy at 408 MHz. 
In addition to the original flux densities at 408 MHz the 
entire classical radio frequency range has been covered for 
most of the sources in the sample, both from already existent 
surveys (6C at 151 MHz, WENSS at 325 MHz, NVSS at 1400 
MHz, GB6 at 4850 MHz) and from dedicated observations, 



mainl y at higher freque ncies (74 M Hz.lMack et al. |2005| ; 2.7 
GHz, Klein et all 120031; 4 85 G Hz, IVigotti et all Il999l and 



10.5 GHzjGregorini et al.ll 19981 ). Based on this radio survey, 
IVigotti et alj (1997:) presented a complete sample of 12EQ 
quasars and this will be used as a comparison sample of 
radio-loud non-BAL quasars. 

In Figure the flux d ensity is plott e d aga inst the red- 
shift but now the sa mple of Vigotti et al, (|l997l ) is presented 
together with that of lBecker et al. ( 200C ) and the one stud- 
ied in this paper. The redshift range covered by the BAL 
QSOs goes from 0.656 to 3.12 4 (Becker's lis t ) and to 3.416 
(RBQ sample). The sources in I Vigotti et al. I i|l997l ) list are, 
on average at smaller redshift (z mea „ = 1.25), but they cover 



1 In practise 124 because the redshift of quasar B3 2329+398 is 
not known. 



3 OBSERVATIONS AND DATA REDUCTION 

Radio-continuum flux densities with full polarisation infor- 
mation were collected for the RBQ sample at different fre- 
quencies in several runs with the Effelsberg 100-m radiotele- 
scope and the Very Large Array (VLA). A summary of the 
different observing runs can be found in Table[2] and the list 
of frequencies used with their characteristics (e.g., angular- 
resolutions) are displayed in Table [3] 



3.1 Effelsberg 100-m telescope 

The single-dish observations at 4.8, 8.4 and 10.5 GHz took 
place in January 2005. Some sources were re-observed in 
July and October 2005 in order to confirm doubtful values 
or improve the signal-to-noise ratio. Observations at 2.65 
GHz were carried out in June 2006. All receivers used in 
these observations are mounted on the secondary focus of 
the 100-m antenna. The recently installed 2.7-GHz hybrid 
feed can observe simultaneously in 9 channels, 8 of them 
10-MHz wide and the ninth covering the full 80-MHz range. 
Since the sources in our sample are too faint to have a high 
signal-to-noise ratio per channel, only the wide-band chan- 
nel centred at 2.65 GHz has been used. The receivers at 
4.8 and 10.5 GHz have multi-feed capabilities with 2 and 
4 horns, respectively, allowing real-time sky subtraction in 
every subscan measurement. 

The observational strategy consisted of cross-scanning 
the source position in azimuth and elevation. These perpen- 
dicular scans were added to gain a factor \/2 in sensitivity 
since the higher resolution FIRST maps (~ 5") showed that 
none of these sources has structure more extended than 1 
arcmin, the higher resolution in our observations (see Table 
[3]). The cross-scan length was chosen to be about four times 
the beam size at each frequency for a correct subtraction 
of linear baselines. The scanning speed ranged from 8" to 
32" /s depending on the frequency, stacking typically 8 sub- 
scans at 2.65 GHz, and some 16 subscans in the rest of the 
frequencies (32 for the faintest sources at 10.5 GHz). This 
translates into on-source integration times between 15-75 sec 
per source and per frequency. Before combining, individual 
subscans were checked to remove those affected by radio fre- 
quency interference, bad weather, or detector instabilities. 
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The calibration sources 3C 286 and 3C 295 were regu- 
larly observed to correct for both time-dependent gain in- 
stabilities and elevation-dependent sensitivity of the antenna 
and t o bring our measu rements to an absolute flux density 
scale (|f3aars et al.lll977T ). The quasar 3C 286 was also used as 
a polarisation calibrator obtaining the polarisation degree m 
and the polaris ation angle y m agreem ent with the values in 
the literature l|Tabara fc Inoug|l980l ) within errors <1 per 
cent and <5 deg, respectively. The unpolarised planetary 
nebula NGC 7027 was also observed to have an estimation 
of the instrumental polarisation at each frequency. 

The measurement of flux densities from the single-dish 
cross-scans was done by fitting Gaussians to the signal from 
all the polarimeter outputs (Stokes I, Q and U) and iden- 
tifying the Gaussian amplitudes with the flux densities Si, 
Sq and Su- For all sources with significant Sq and Sv con- 
tributions, the polarised flux density Sp, the degree of linear 
polarisation m and the polaris ation angle y were computed, 
using the expressions given in lKlein et alj (|2003T ). 



3.2 Very Large Array 

VLA data were taken during February/March 2006 in the 
most extended A configuration using the receivers corre- 
sponding to the X, U, K and Q bands, i.e., at 8.4, 15, 22 
and 43 GHz, respectively. Integration times were chosen on 
the basis of the predicted flux densities at these frequencies, 
after extrapolating from the Effelsberg frequencies with the 
adequate spectral indices. Those sources with indications of 
measurable polarised emission were observed more deeply 
and those expected to be too faint for this purpose were ob- 
served to have at least a good signal to noise ratio in Si. 
Integration times range from ~1 minute at 8.4 GHz up to 
~1 hour at 22 or 43 GHz for some sources. 

At the two highest frequencies the fast-switch method 
was used, quickly switching between target and calibrator 
in order to have phase stability, taking care of the rapid 
atmospheric fluctuations. The duration of the cycles target- 
calibrator was of ~ 150 to 200 seconds in the two highest 
frequencies. 

3C 286 was the primary flux density calibrator and 
3C 48 was used to test the goodness of the flux scale so- 
lution. Several phase calibrators chosen from the VLA cal- 
ibrator manua0 were observed during the run. Most were 
selected at about 2-5 deg from their target sources, and in 
all cases within 10 deg. This is especially important to avoid 
the loss of coherence at the h ighest frequencies. Again the 
flux density scale is the one of|B aars et all l|l977h except at 
22 and 43 GHz where the Baars expressions are not valid. 
At these frequencies the adopted scale is based on emission 
models and observations of planet Mars, and is supposed to 
be accurate at a level of 5-10 per cent. 

The 31DEC06 version of the AIPS package was used to 
flag, calibrate, clean and image the VLA data. The standard 
recipe was followed in all frequencies. However, at 22 and 43 
GHz opacity and gain curve corrections were applied as a 
previous step to the standard calibration, and the calibra- 
tion table was gridded in intervals of 3 seconds to get rid of 



2 Available at |http://www. via. nrao.edu/ast ro/calib/manual 
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Table 3. Observing frequencies and beam sizes. 



Telescope 


Frequency 


Bandwidth 


%PBW 




(GHz) 


(MHz) 


(arcscc) 


Effelsberg 


2.65 


350 


265 


Effelsberg 


4.85 


500 


145 


Effelsberg 


8.35 


1200 


80 


Effelsberg 


10.5 


300 


65 


VLA(A) 


8.45 


700 


0.24 


VLA(A) 


14.5 


700 


0.14 


VLA(A) 


22.5 


2000 


0.08 


VLA(A) 


43.5 


10000 


0.05 



rapid atmospheric fluctuations. The appropriate clean com- 
ponent models were used for 3C 286 and 3C 48 at these two 
frequencies. The peak flux densities were measured by fit- 
ting 2-D Gaussians to the source profiles. The integrated flux 
densities were computed using the BLSUM task, adding the 
signal in small boxes containing each source. The same boxes 
were used to measure the Stokes parameters Sq and Su- The 
local noise of the maps was estimated from a circular region 
around each individual source. Peak and integrated flux den- 
sities have essentially the same values because most sources 
are unresolved, as will be discussed below. 

The polarisation calibration of the VLA data was de- 
termined using the AIPS tasks PCAL and RLDIF. The first 
solves for the feed parameters in each antenna, using re- 
peated observations of a strong unresolved source (in our 
case several snapshots of 3C 286 were taken) covering a range 
of different parallactic angles. RLDIF determines and sub- 
tracts any phase difference between the right and left po- 
larisation systems. The residual instrumental polarisation 
at the VLA frequencies was estimated measuring the linear 
polarisation of 3C 84 a f ter h aving applied the polarisation 
calibration. lAller et alj (|2003T ) report that this source shows 
a percentage of linear polarisation of 0.12 ± 0.01 at 14.5 
GHz, and for our purposes it can be considered unpolarised. 

3.3 Error determination 

We consi der three ma i n cont ributions to the flux density er- 
ror as in iKlein et all l|2003h . These are (i) the calibration 
error AS ca i (in percentage) which is estimated by the dis- 
persion of the different observations of the flux density cal- 
ibrators; (ii) the error introduced by noise, AS n ,i, which 
is estimated from the local noise around the source; and 
(iii) the confusion error AS C onf,i due to background sources 
within the beam area. This last term can be neglected in 
interferometric data where the synthesised beam has small 
dimensions. 

Thus the expressions for the total uncertainty of Stokes 
parameters are Equations [1] and [2] for Effelsberg and the 
VLA, respectively: 

ASi = y/(Si ■ AS cal ) 2 + ASl t + AS 2 conU i = I,Q,U(l) 



AS t = J {Si ■ AS cal ) 2 + Sl ti ■ i = I, Q, u (2) 

where A src is the area of the box used to measure the 
source flux density, and Abeam is the beam area. The ex- 
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Table 4. Values used to calculate the errors of the 
Stokes parameters and obtained instrumental polari- 
sations. 



Telescope Frequency A5 con/ AS cai Pinstr 





^Lrrlz ) 


(1) 


(,q,uj 


(per 


cent) 


Effelsberg 


2.65 


1.5 


0.5 


0.8 


0.7 


Effelsberg 


4.85 


0.45 


0.15 


0.9 


0.7 


Effelsberg 


8.35 


0.23 




1.0 


0.6 


Effelsberg 


10.5 


0.08 




1.9 


0.8 


VLA 


8.45 






1.9 


<0.1 


VLA 


15.0 






2.9 


0.6 


VLA 


22.5 






1.8 


0.7 


VLA 


43.5 






5.3 


2.1 



press ions for the uncert ainties of Sp, m and x can De taken 
from iKlein et al l (|2003h . 

Table [4] shows the coefficients used for the calculation 
of errors AS;. The confusio n limits at 2 . 7, 4.8 and 10.5 GHz 
have been extracted from [Klein et al.l ll200St ) and that at 
8.35 GHz (not available in lKlein et ail 120031 ) was obtained 
from our own observations inspecting the behaviour of the 
signal rms at long integration times. This value is consistent 
with the interpolation over frequency from the two adjacent 
bands (at 4.85 and 10.5 GHz). The last column in Table g] 
lists the instrumental polarisation for the different frequen- 
cies and telescopes. 



4 RESULTS 

4.1 Radio flux densities 

Tables [5] and [6] show the measured flux densities and errors 
(in mjy) between 74 MHz and 43 GHz for the RBQ sample. 
At frequencies below 2.6 GHz there are available flux densi- 
ties from the literature and these are listed in Table [5] The 
flux densities at 325 MHz were obt ained from the Wester - 
bork Northern Sky Survey f WENSS; iRengelink et af1ll997h . 
In addi tion, 3-a upper limit s and a detection from the Texas 
survey (iDouglas et al.lll996T ) at 365 MHz are given for those 
sources not covered by the WENSS survey. At 74 MHz none 
of the sources is detected in the VLA Low-Frequency Sky 
Survey (VLSS; ICohen et al.ll2007h . but 3-<t upper limits have 
been computed measuring the local noise around the source 
position in the available i mages. Flux densitie s at 1.4 GHz 
from FIRST and NVSS (|Condon et all Il998l ) can be also 
found in Table [5] 

The BAL QSO 1624+37 was observed in this campaign 
at 2.65 GHz with Effelsberg and at 15 and 43 GHz with 
the VLA. All the flux densities in th e remaining frequencies 
were presented bv lBenn et all l|2005h . 

4.2 Morphologies 

The FIRST maps at 1.4 GHz show a compact morphology 
for all 15 BAL QSOs in the RBQ sample, at a resolution 
of ~5". To quantify this we have computed the morpholog- 
ical parameter O, based on the FIRST integrated and peak 
flux densities, O = log(Si n t/S pea k), and for all 15 sources 
<0.03. All sources, with the exception of 1053—00, appear 



Table 5. RBQ sample of 15 radio-loud BAL QSOs: Flux den- 
sities (in mjy), errors or 3-cr upper limits from 74 MHz to 1.4 
GHz obtained from the literature. 



TD cVLSS ct cFIRST cNVSS 

1JJ °0.074 ^0.365 °1.4 "1.4 



0039-00 


<475 


<80 


21.2 


19.5 ± 0.7 


0135-02 


<350 


<80 


22.4 


22.8 ± 1.1 


0256-01 


<385 


<80 


27.6 


22.3 ± 0.8 


0728+40 


<300 


15 ± 3.3 1 " 


17.0 


17.2 ± 0.6 


0837+36 


<290 


184 ± 3.4t 


25.5 




0957+23 


<260 


<80 


136.1 


136.9 ± 4.1 


1053-00 


<500 


<80 


24.7 


22.7 ± 1.1 


1159+01 


<320 


887 ± 30 


268.5 


275.6 ± 8.3 


1213+01 


<340 


<80 


21.5 


27.5 ± 7.9 


1228-01 


<340 


<80 


29.4 


29.1 ± 1.0 


1312+23 


<300 


<80 


43.3 


46.5 ± 1.4 


1413+42 


<330 


16 ± 3.2t 


17.8 


16.8 ± 0.6 


1603+30 


<225 


33 ± 4.4t 


53.7 


54.1 ± 1.7 


1624+37 


<355 


59 ± 4.0 1 " 


56.1 


55.6 ± 1.7 


1625+48 


<290 


26 ± 3.7t 


25.3 


26.0 ± 0.9 



The symbol f indicate data from WENSS at 325 MHz, the 
rest are from the Texas Survey at 365 MHz. 



essentially point-like at this resolution. The FIRST map of 
1053—00 is shown in Figure[3]where a point-like core and two 
faint components (5- and 8-cr detections) can be seen. The 
faint components seem to be real and this source has been 
includ ed in the catalogue of double-lobed radio quasars from 
SDSS (|de Vries et all 120061 '). This is thus an additional ex- 
ample of the rare class of BAL QSOs showing double-lobed 
morphology, and it is not inclu ded in the list of 8 S DSS BAL 
QSOs of this type compiled bv lGregg et ail i|2006f ). 

Our VLA maps at higher frequencies confirm the com- 
pactness of the 15 BAL QSOs, having all of them point-like 
structure at 8.4 and 15.0 GHz. The faint components of 
1053—00 cannot be seen at these frequencies, probably due 
to a steep spectral index in these two regions. 

The combination of the 43-GHz band with the A con- 
figuration provides the highest resolution available with the 
VLA {Oh p bw ~ 50 mas). However, as can be seen in Table 
[Sjmost sources are very faint at 43 GHz, i.e., close to the 1 
mjy level. For those the self-calibration process is not ap- 
propriate and the way to recover the flux densities was to 
assume an initial point-source model before the first map- 
ping. This procedure is to some extent well justified since 
at lower frequencies the morphology is still point-like, al- 
though the precise astrometric and morphological informa- 
tion about the sources at 43 GHz is lost. 

Self-calibration at high frequency was only possible for 
the brightest sources: 0957+23, 1159+01 and 1312+23. The 
first two sources appear point-like in the 43-GHz maps. Fig- 
ure [4] shows the naturally weighted maps of 1312+23 at 43 
and 22 GHz, revealing an extension towards the north com- 
ing from the central unresolved component. Although the 
feature, possibly a jet, is detected at both frequencies, its 
significance is higher at 43 GHz. 

The small feature in the 43-GHz map at ~120 mas south 
from the central source is probably not real because it is lo- 
cated in a region affected by a stripe. It must be noticed that 
at this high frequency this flux density level is probably close 
to the limit where signal coherence allows accurate mapping. 
For the remaining 12 sources the mapping at 43 GHz was not 
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Table 6. RBQ sample of 15 radio-loud BAL QSOs: Flux densities (in mjy) from 2.6 to 43 GHz obtained in this work. Superscripts on 
the errors indicate the run number of that ob servation, accordin g to the key in Table [2] The exception is 1624+37 (with superscript 0), for 
which several flux densities were presented in lBenn et al ] l|2005ri 



ID 




52.6 




S4.8 




Sg.3 




Ss.4 




SlO.5 




S15 




S22 


S43 


0039-00 




— 




19.3 


+ 


0.6 1 


12.0 


± 


0.4 1 


12.3 


+ 


0.1 5 


9.7 ± 


Q 4l 


6 


2 


+ 


l.l 5 


2.2 ± 1.2 5 


1.1 ± 0.1 5 


0135-02 




— 




27.2 




0.6 1 


31.9 




0.8 1 


30.8 


4- 


0.4 5 


29.5 + 


0.7 3 


22 


7 




0.8 s 


14.7 ± 0.2 5 


2.3 ± 0.2 5 


0256-01 




— 




12.0 


+ 


0.5 1 


7.2 


+ 


0.3 1 


7.2 


+ 


0.4 5 


5.9 ± 


0.4 1 


1 


3 


+ 


1.5 s 


2.2 ± 0.3 5 


1.4 ± 0.1 5 


0728+40 


6 


7 ± 


1.5 4 


6.4 


+ 


0.7 1 


3.1 


+ 


0.2 1 


2.5 


+ 


0.3 5 


1.6 ± 


31 


1 


9 


+ 


0.3 5 


0.83 ± 0.06 5 


— 


0837+36 


50 


4 ± 


1.6 4 


30.2 


± 


2.3 2 


12.1 


± 


0.8 1 


12.5 


± 


0.1 6 


13.0 ± 


0.4 2 


3 


3 


± 


1.3 6 


2.2 ± 0.5 6 


1.3 ± 0.2 6 


0957+23 


93 


7 ± 


2.4 4 


78.7 


± 


l.l 1 


48.5 


± 


1.2 1 


53.2 


± 


0.5 5 


39.8 ± 


1.2 1 


31 


4 


± 


0.7 5 


22.4 ± 0.3 5 


9.8 ± 0.9 5 


1053-00 


13 


6 ± 


1.6 4 


15.7 


± 


0.6 1 


11.8 


± 


0.4 1 


13.1 


± 


0.1 5 


12.3 ± 


0.6 1 


8 


2 


± 


1.6 5 


5.9 ± 0.7 5 


2.3 ± 0.7 5 


1159+01 


133 


9 ± 


2.0 4 


137.8 


± 


1.7 1 


158.0 


± 


2.0 2 


160.8 


± 


1.2 5 


150.6 ± 


3.7 2 


123 


8 


± 


1.6 5 


105.5 ± l.l 5 


74.6 ± 2.1 5 


1213+01 


23 


± 


1.7 4 


15.2 


± 


0.5 1 


9.7 


± 


0.3 1 


11.6 


± 


0.1 5 


8.2 ± 


0.5 1 


5 


6 


± 


1.4 s 


3.6 ± 0.5 5 


1.4 ± 0.2 5 


1228-01 


19 


4 ± 


1.6 4 


18.6 


± 


0.5 1 


15.5 


± 


0.4 1 


16.9 


± 


0.2 s 


13.5 ± 


0.7 1 


11 


5 


± 


1.6 5 


8.1 ± 0.2 5 


2.5 ± 0.7 5 


1312+23 


28 


9 ± 


1.6 4 


25.7 


± 


0.6 1 


19.8 


± 


0.5 1 


19.4 


± 


0.2 s 


16.1 ± 


0.6 1 


10 


5 


± 


0.3 5 


7.7 ± 0.1 6 


5.4 ± 0.4 6 


1413+42 


7 


6 ± 


2.9 4 


8.8 


± 


0.7 2 


13.4 


± 


0.3 2 


13.4 


± 


0.1 6 


12.7 ± 


0.4 2 


9 


5 


± 


1.5 6 


7.0 ± 0.2 6 


2.2 ± 0.4 6 


1603+30 


22 


8 ± 


1.7 4 


26.1 


± 


0.7 1 


19.1 


± 


0.5 1 


22.1 


± 


0.3 5 


16.8 ± 


0.6 1 


12 


5 


± 


1.7 s 


8.3 ± 0.1 5 


1.9 ± 0.5 5 


1624+37 


33 


5 ± 


1.6 4 


23.3 


± 


1.1° 








15.0 


± 


0.1° 


10.5 ± 


0.8° 


9 


6 


± 


0.6 s 


5.41 ± 0.02° 


2.1 ± 0.3 5 


1625+48 


17 


6 ± 


1.5* 


9.4 


± 


0.7 2 


8.0 


± 


0.2 2 


7.0 


± 


0.2 5 


6.3 ± 


0.2 2 


5 


8 


± 


1.3 s 


1.9 ± 0.3 5 


1.7 ± 0.2 5 



1053-00 Total intensity 1.4 GHz 

i i 



Table 7. Angular dimensions for the extended sources. 








10 53 55 54 53 52 51 

RIGHT ASCENSION (J2000) 

Levs = 1 .6E-04 * (-3, 3, 4, 5, 8, 1 6, 32, 64, 1 28) mjy/beam 



Figure 3. Map of BAL QSO 1053-00 from FIRST at 1.4 GHz. 
The synthesized beam-size is shown in the lower-left corner. 



possible and an upper limit on their sizes was obtained using 
the observations at 22 GHz (Qhpbw ~ 80 mas), where they 
are bright enough for mapping. For those, the task IMFIT 
was used to fit an elliptical Gaussian component and a zero 
level within a small box containing the radio source, yielding 
the major and minor axes of the fitted Gaussian, position 
angle and the formal beam-deconvolved dimensions. The un- 
certainty in_the_dimensions was estimated following Equa- 
tion 1 o f lWhite et ail (|l997l V Eight of the twelve sources are 
clearly more extended than the beam size and their fitted 
parameters are given in Table[7] Also in this table are shown 
the dimensions of 1312+23 and 1624+37, as measured from 
the map in Figure|4]and from VLBA observations to be pub- 
lished in a subsequent paper, respectively. The remaining 



Source 


Freq. 
(GHz) 


(mas) 


(mas) 


PA 

(deg) 


n(dec) 
maj 

(mas) 


0135-02 


22 


112.4 


78.7 


-12 


11 + 4 


0728+40 


22 


85.5 


73.0 


26 


26 + 6 


0837+36 


22 


94.0 


85.9 


34 


19 ± 10 


1213+01 


22 


111.7 


75.5 


-19 


27 ± 11 


1228-01 


22 


95.3 


76.3 


-30 


19 + 4 


1312+23 


43 


54.5 


41.8 


61 


69 ± 12t 


1413+42 


22 


92.0 


69.5 


-60 


9 + 3 


1603+30 


22 


95.8 


71.0 


-71 


15 + 3 


1624+37 


4.8 








36 ± 5tt 


1625+48 


22 


157.0 


71.8 


76 


63 ± 20 



t Not deconvolved size, but measured from the extended 
structure in Figure|4] ft Measured from VLBA map at 4.8 GHz 
to be presented in Montenegro-Montes et al. (in prep.) 



three sources (0039-00, 0256-01, 1053-00) plus 0957+23 
and 1159+01 can be considered strictly unresolved at high 
frequencies. 

VLBI obs ervations of 095 7 +23 a nd 1312+23 have been 
presented by Ijiang fc Wanel l|2003h using the European 
VLBI Network at 1.6 GHz. With a restoring beam of only 
18.5 x 6.86 mas, 0957+23 appears as a single point-like com- 
ponent, while 1312+23 is resolved. This last source shows a 
symmetric morphology with a central core and two compo- 
nents in the northern and southern directions with a total 
extension of ~150 mas (i.e., about 1 kpc). Although an ex- 
tended morphology is also present in our 22-GHz and 43- 
GHz maps (Fig. [4}, only the northern part is detected. 

Summarising, our VLA maps show very compact mor- 
phologies for all 15 sources at all frequencies, at a maximum 
resolution of 80 mas. The exceptions are 1312+23 which 
shows some elongation at 22 and 43 GHz, and 1053—00 
which is extended at 1.4 GHz showing two faint lobes not 
detected at higher frequencies. 
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1312+23 Total Intensity 22460 MHz 

i i i i r 




Table 8. Sources showing significant flux density variability. At 
1.4 GHz, Si co mes from NVSS and S 2 from FIRST. At 8.4 GHz, 
Si comes from lBecker et ah | (|200(]| 1. S 2 from this work (Table© 



Source 


Freq 


Si 


S 2 


Var^s 


0~Var 




(GHz) 


(mjy) 


(mjy) 






0256-01 


1.4 


22.3+0.8 D 


27.5 B 


0.23 


3.2 


1213+01 


1.4 


27.5+0.9 D 


22.9 B 


-0.20 


3.2 


1312+23 


8.4 


12.6 A 


19.4+0.2 A 


0.54 


10.3 


1413+42 


8.4 


11.3 D 


13.4+0.1 A 


0.19 


3.7 


1603+30 


8.4 


18.1 A 


22.1+0.3 A 


0.22 


4.2 



58.2 - 



13 12 13.60 13.59 13.58 13.57 13.56 13.55 
RIGHT ASCENSION (J2000) 

Levs = 6.000E-03 (-3, 3, 4, 5, 8, 16, 32, 64) mjy/beam 
1312+23 Total Intensity 43340 MHz 

23 19 58.85 

58.80 

58.75 

58.70 

i 58.65 

' 58.60 

| 58.55 
I 
I 

58.50 
58.45 
58.40 
58.35 




o 



13 12 13.595 13.590 13.585 13.580 13.575 13.570 13.565 13.560 13.555 
RIGHT ASCENSION (J2000) 

Levs = 1.250E-01 (-3, 3, 4, 5, 8, 16, 32, 64) mjy/beam 



Figure 4. VLA map of BAL QSO 1312+23 at 22 GHz (top) and 
at 43 GHz (bottom). The synthesized beam-size is shown in the 
lower-left corner of each map. 



4.3 Variability 

We have checked for variability at 1.4 GHz comparing the 
flux densities from FIRST and NVSS epochs in those 14 
sources with NVSS measurement. The same has been done 
at 8.4 GHz for tho se 5 sources in common with the list of 
iBecker et al.l ((2000). They observed these sources with the 
VLA in A and D configuration (resolutions of 0.8 and 9 
arcsec respectively) and all 5 sources presented point-like 
structure at both resolutions. As mentioned before this is 
again confirmed by our observations at 8.4 GHz in A con- 
figuration. This means that flux densities of both measure- 
ments can be directly compared. The same is true at 1.4 
GHz because all sources except 1053—00 are point-like at 
the resolution of FIRST (B configuration) and the resolu- 



tion of NVSS is poorer (VLA D and CnB configurations). 
Although 1053—00 is resolved, it has a compact nucleus and 
the two extensions are not expected to contribute much to 
the total flux density (Section 14. 2[) . 

We measure the flux densi ty variations with the p aram- 
eter VarAS, as denned e.g. by Torniainen et al l <|2005h : 



VarAS 



Sri 



Sri 



Sri 



(3) 



As an estimate of the significa nce of the source v ariability 
the a var parameter defined by IZhou et al.1 (120061 ) is used, 
adopting the integrated flux densities Si and S2 at the two 
epochs: 



0~Var = 



\S 2 -Si\ 



(4) 



where 02 and tj\ are the uncertainties in the integrated flux 
densities. Table [S] shows flux densities, Var^s and ovar at 
1.4 and 8.4 GHz, for those sources exhib iting avar > 3. For 
the flux densities of IBecker et all (|2000l ). whose errors are 
not given in their work, a 5 per cent error is assumed. 

A low percentage of the sample (only 2 out of 14 
sources) shows variability at 1.4 GHz, both sources vary- 
ing about 20 per cent and with a significance just above 
3 standard deviations. At 8.4 GHz, 3 out of 5 sources are 
variable at this level and 1312+23 constitutes the most ex- 
treme case, showing variations of about 50 per cent with 
significance above 10. 

Flux density variations have recently been used to con- 
strain the orient ation of several BAL Q SOs. IZhou et al.l 
l|2006l ) and later iGhosh fc Punslvl (|2007l ) found a few ex- 
amples of BAL QSOs with variability as high as 40 per cent 
at 1.4 GHz between the NVSS and FIRST, i.e., in periods of 
about 1-5 years. These variations were high enough to make 
the brightness temperature Tb exceed the 10 12 K limit asso- 
ciated to the Compton catastrophe. This was interpreted as 
a sign of beaming and the jet was supposed to be oriented 
at a few degrees from the line of sight. 

The same procedure can be applied to 0256—01 and 
1213+01, accept ing the 3-tr signif i cance as real. Following 
the formalism in lGhosh fc Punslvl l|2007l ) and knowing that 
the time differences between the NVSS and FIRST observa- 
tions are 2.1 and 3.4 years, respectively, the observed vari- 
ations imply brightness temperatures of 21=3. 1T0 13 K for 
0256-01 and 1.3-10 13 K for 1213+01. These translate into 
upper limits for the jet orientation angle with respect to the 
line of sight of ~20 and 25 deg, respectively. One concern 
about this approach is that it is based on variability at a sin- 
gle frequency and from only two measurements. Moreover, 



© 2008 RAS, MNRAS 000,rflfl8l 



Radio spectra and polarisation of BAL QSOs 9 



the flux densities we obtained for 0256-01 at 8.4 GHz, with 
Effelsberg and almost one year later with the VLA, show a 
good match, i.e., 7.2±0.3 and 7.2±0.4 mjy respectively, in 
contradiction with the expected stronger variability at high 
frequencies. 

4.4 Shape of the radio spectra 

4-4-1 Sample of radio BAL QSOs 

Figure [5] show the spectra of the 15 objects. Many of them 
have a relatively steep spectrum at high frequency. About 75 
per cent of the sample show flattening of the spectrum at low 
frequencies (typically below 1-5 GHz) which could indicate 
synchrotron self-absorption. However, in some cases some 
variability cannot be excluded, since the 1.4-GHz measure- 
ments were taken some years before the others. About 1/3 of 
the sample show, in addition, enhanced emission at MHz fre- 
quencies which could be indications of a second component 
emitting mostly in this range. Again, variability cannot be 
ruled out, but any possible time variability is expected to be 
less pronounced at lower frequencies. The source 1159+01 is 
the one which best fits in this category with both, a depres- 
sion of the spectrum towards low frequency and a contribu- 
tion from a second component i n its spectrum. 

Following the approach by iDallacasa et ail (|2000T ). we 
have fitted each spectrum with the analytical function: 

logSu =a - \Jb 2 + (c log(u) - d) 2 (5) 

where a, b,c,d are numeri c al con stants. This function, pre- 
sented in IDallacasa et all (|200Ch was used by these authors 
to obtain the frequency turnover, i^peak, of sources with con- 
vex spectra peaking above a few GHz, referred to in that 
paper as High Frequency Peakers. The value of !/ pea k ap- 
proximately marks the frequency region where the source 
starts to be optically thick to synchrotron radiation. 

These fits have been done excluding some points in or- 
der to obtain a reduced x 2 below 0.05. This was possible for 
most of the sources which present a smooth radio spectrum, 
specially at high frequency, but not for 0728+40 (x 2 =0.08) 
which presents a large scatter. At low frequencies, as we 
noted before, there are sources for which a possible second 
component arises. Again in these cases we decided whether 
to include or not in the fit these (low frequency) points try- 
ing to maintain the reduced x 2 within the mentioned limit. 
An exception was done with 0957+23 because even if an ac- 
ceptable \ 2 was obtained when including the 1.4 GHz flux 
densities, the resulting fit was not compatible with the 365 
MHz upper limit. 

In our sample we find 9 sources for which the spectrum 
peaks at a frequency higher than 1 GHz in the observer's 
frame and there are 4 displaying a peak in the range 300-700 
MHz. For the remaining three objects (0256-01, 0957+23 
and 1312+23) the spectrum is slightly convex but not show- 
ing an obvious peak. However, power-law fits would not be 
compatible with the Texas upper limits and then extrapo- 
lated peaks compatible with both the spectrum curvature 
and these upper limits have been determined. The peak fre- 
quencies have to be multiplied by (1+z) to obtain intrinsic 
peak-frequencies. 

In order to understand the global trends on the spectral 
shape of BAL QSOs, we compare in Figure [6] all the spec- 



tra bringing them first to the rest-frame in order to avoid 
the effect of the redshift. The rest-frame normalisation fre- 
quency was chosen to be 25 GHz. This is a compromise to 
choose a region representative of pure synchrotron emission, 
where neither synchrotron self-absorption nor synchrotron 
losses might be present in the majority of sources. This fre- 
quency also makes the dispersion of the normalised spectra 
relatively small compared to other normalisation frequen- 
cies. 

We see in Figure [S] a variety of spectral slopes around 
25 GHz. There is one object in the sample which differ more 
strongly from the main trend, This is 1413+42 which has a 
spectrum with two components separated at 4 GHz (see Fig- 
ure [5J, and for this reason it shows in Figure [6] an inverted 
spectrum around the normalisation frequency. Apart from 
this object, and assuming that variability effects are not 
very important, Figure [6] reflects the rest-frame spectrum of 
a typical BAL QSO in the RBQ sample. It is in general con- 
vex, quite flat below 10 GHz, and becoming steeper around 
25 GHz. About half of the objects become even steeper at 
frequencies higher than 50 GHz. 

The spectral shape can be quantified by means of the 
spectral indices. Table [5] provides the turnover frequencies 
and representative spectral indices at different frequencies, 
both in the observer's and in the rest-frame. The spectral 
indices ao!365> a 4% an d a| 3 4 describe the low-, medium- and 
high-frequency r egions. Moreover, th e index af;| is shown for 
comparison with iBecker et all (|2000h . It is worth noting that 
the spectral indices af;| shown in Tableware not based on 
simultaneous observations, but as we previously have shown, 
variability might not be very important at this frequency, 
and in any case it only affects a small number of sources. In 
the rest-frame the indices al 5 o, a 2 |, ali and 050° are also 
given in Tabled 

Looking at the rest-frame spectral indices of BAL QSOs 
and Figure [5] we see on average a flattening below 10 GHz, 
which is probably due to synchrotron self-absorption, and 
an overall convex shape. The low-frequency flattening is 
reflected in the median a§ 5 =— 0.32, while the median 0^2 
decreases to —0.57. The steepening at high frequency can 
be seen from the median a^s dropping to —0.92 and fi- 
nally becoming steeper at higher frequencies with a median 
«|8°=-1.24. 

The radio spectral index is known to be a statistical 
indic a tor of the orientatio n of radio sources (|Orr fc Browne! 
1982). Becker et al. (2000) found respectively two thirds and 
one third of the radio-loud BAL QSOs in their list with steep 
and flat spectral index. However, it is worth to note that they 
define the limit between flat and steep to be a=— 0.5, and 
about one third of the sample has spectal index close to this 
value —0.6 a ^ —0.4. They suggest that this result is not 
consistent with BAL QSOs being oriented along a particular 
direction with respect to the line of sigh t to them. We find in 
our sample, which overlaps with that of lBecker et alj (|2000h 
in 5 objects, 60 per cent (9/15) of flat and 40 per cent of steep 
(6/15) spectrum sources, with 27 per cent of objects (4/15) 
being in the range —0.6 ^ a ^ —0.4. These percentages are 
roughly consistent within the errors, which are large because 
the samples are relatively small. 
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Figure 5. Radio spectra of the 15 BAL QSOs in the RBQ sample. Errors are shown if larger than the symbol size. Triangles mean 3-cr 
upper limits. Solid lines are fits to the analytical function give in Section l4,4.1l The peak frequencies given at the top- right of each pannel 
were obtained from the fit. Some sources like 1159+01, 1413+42 and 1603+30 seem to have a second component at low frequencies. 
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Table 9. Turnover frequencies (in GHz) and spectral indices for the RBQ sample in the observer's frame 
(columns 2-6) and the rest frame (columns 7-11). The spectral indices in the observer's frame have been 
computed from each pair of observed flux densities while the rest-frame spectral indices are based on 
interpolation over the whole spectrum. Median values of the spectral indices are given in the last row. 



ID 

(1) 


Speak 

(2) 


Observer's frame 
"0.365 "1.4 "4.8 

(3) (4) (5) 


"! 3 4 

(6) 


Speak 

(7) 


"| 5 o 
(8) 


Rest frame 

«?§ 
(9) 


< 
(10) 


(11) 


0039-00 


2.60 




-0.30 


— 1.01 


— 1.48 


8.45 


-0.32 


-0.54 


—1.21 


-1.88 


0135-02 


8.15 


— 


0.18 


-0.16 


-1.59 


22.90 


+0.17 


+0.18 


-0.70 


—2.26 


0256-01 


1.60 




-0.75 


-0.91 


-1.01 


5.50 


-0.74 


—0.80 


—0.90 


— 1.28 


0728+40 


0.70 


+0.08 


-1.07 


-1.08 


— 


1.15 


-0.88 


-0.75 






0837+36 


2.70 


-1.35 


-0.40 


-1.96 


-1.39 


12.00 


—0.09 


-1.00 


-1.30 


—2.16 


0957+23 


1.40 




-0.52 


-0.81 


-1.04 


4.25 


-0.51 


-0.60 


-0.91 


-1.09 


1053-00 


4.40 




-0.35 


-0.58 


-1.07 


11.15 


-0.14 


-0.30 


-0.94 


-1.31 


1159+01 


6.30 


-0.81 


-0.29 


-0.09 


—0.47 


18.90 


—0.09 


+0.22 


—0.44 


—0.48 


1213+01 


1.30 




-0.35 


-0.88 


-1.30 


4.90 


-0.35 


-0.61 


-1.02 


-1.20 


1228-01 


4.15 




-0.31 


-0.43 


-1.17 


15.10 


-0.29 


-0.12 


-0.55 


-1.12 


1312+23 


0.35 




-0.45 


-0.79 


-0.79 


0.85 


-0.43 


-0.58 


-1.01 


-0.57 


1413+42 


9.20 


+0.07 


-0.16 


+0.07 


-1.11 


35.05 


-0.23 


+0.45 


-0.06 


-1.04 


1603+30 


4.50 


+0.33 


-0.50 


-0.65 


-1.51 


13.50 


-0.29 


-0.16 


-0.98 


-1.73 


1624+37 


0.55 


-0.03 


-0.74 


-0.79 


-1.21 


2.40 


-0.71 


-0.64 


-0.99 


-1.00 


1625+48 


0.65 


-0.02 


-0.72 


-0.43 


-0.87 


2.45 


-0.76 


-0.81 


-0.41 


-1.69 


Median: 




-0.02 


-0.40 


-0.79 


-1.14 




-0.32 


-0.57 


-0.92 


-1.24 



4-4-2 Comparison with non-BAL QSOs 

To further investigate the question of the orientation of BAL 
QSOs we want to compare in this paper the spectral index 
distribution of quasars in a pure BAL QSO sample with that 
of a sample of non-BAL QSOs with the idea of interpret- 
ing possible differences with different orientations of both 
samples. In this comparison we will use the spectral range 
1.4-8.4 GHz in the observer's frame and 6.0-25 GHz in the 

rest frame. 

The list of non-BAL QSOs will be that of IVigotti et all 
l| 19971 ) presented in Section [2] but in all the following discus- 
sion we will restrict it to objects with z>0.5. This restriction 
only reduces the total number of objects from 124 to 119. 
We have computed the spectral indices in the same way as 
in our sample, considering pairs of observed flux densities in 
the observer's frame and interpolating over the whole spec- 
tra in the rest frame. The spectral range covered by the 
non-BAL QSO sample goes from 74 MHz up to 10.5 GHz, 
and only those sources with z>1.5 will reach 25 GHz in the 
rest frame. 

As a BAL QSO sample w e will add to our RBQ sample 
the list of (jBecker et all 12000). They only provide flux den- 
sities in two frequencies, i.e., 1.4 and 8.4 GHz. Thus, com- 
bining these we end up with a combined sample of 38 BAL 
QSOs. These spectral indices can be found i n Table [9] for 
objec ts in the RBQ sample, and in Table 2 of dBecker et al.l 
120001 ) for the remaining objects. 

The Vigotti et al. sample could in principle contain both 
BAL QSOs and non-BAL QSOs . However, a low con tamina- 
tion of BAL QSOs is expected. iBecker et al.l (|200lh showed 
that in the radio-powerful regime the fraction of BAL QSOs 
is small. They u sed the defini t ion of radio-loudness param- 
eter R* given bv lStocke et~ai1 (|l992T ), i.e., the ratio of radio 
power at 5 GHz to optical luminosity at 2500 A in the rest 
frame. In the regime of QSOs with log(R*) >2 they esti- 
mated a fraction of 4 per cent (<1.5 per cent) of HiBAL 
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Figure 6. Rest-frame spectra of 15 BAL QSOs in the 
RBQ sample. They have been normalised to v r est = 25.0 GHz 
(log(iv est )=1.398) 

(LoBAL) QSOs from the FBQS, in the interval of redshift 
[1.5,4.0] ([0.5,1.5]) where the Civ (Mgn) absorbing features 
would be covered in optical spectra. 

We have calculated log(R*) for all quasars in Vigotti's 
list. For this purpose, the radio luminosity at 5 GHz has 
been computed on the basis of the observed radio spectra, 
applying the appropriate K-correction interpolating between 
each pair of observed frequencies with the adequate spec- 
tral index. The optical luminosities at 2500 A were com- 
puted using the PO SS-I APM E magnitudes published in 
IVigotti et al](|l997n , but corrected from Galactic extinction 
following the maps of lSchlegel et all l| 19981 ) and assuming an 
optical spectral i ndex of a ovt = — 1 ( S v oc v a ), which is the 
value adopted bv lBecker etail l|200ll ). 
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Fi gure 7. Histogram of radio-loudness parameter R* as defi ned 
bv lStocke et al.l lll992h for the sample of lVieotti et al.l i 19971) . 



Figure 8. Spectrum of SDSS J080016.09+402955.6, showing a 
relatively narrow absorption in C IV and Si IV. 



In Figure [7] the histogram shows how log(R*) ranges 
from 1.9 to 4.6. To estimate the fraction of BAL QSOs fol- 
lowing the mentioned proportions, we split this sample in 
two groups, one with 0.5<z^l.5 and the other with z>1.5, 
comprising 80 and 44 quasars respectively. About 1-2 Hi- 
BAL QSOs are expected in the high-z subsample and about 
1 BAL QSO in the low-z subsample. 

We made use of the NASA Extragalactic Database 
(NED) and the SDSS-DR 6 database to l ook f or the opti- 
cal spectra of QSOs from IVigotti et alj l| 19971 ) and check 
whether they are BAL QSOs. From 119 objects with z<0.5, 
55 were found in the SDSS-DR6 database. These were in- 
spected visually to look for the presence of BAL features. All 
of them are non-BAL QSOs and only one of them, (SDSS 
J080016.09+402955.6, associated to the radio source B3 
0756+406) shows some broad absorption component blue- 
wards of the Civ emission line, also present in Siiv. How- 
ever, this absorption has a width of ~ 1000 km s _1 and 
it is located within 3000 km s _1 from the line emission. It 
is therefore not consistent with the definition of BAL QSO 
but it could be better classified as a mini-BAL or a Narrow 
Absorption Line (NAL) system (see Figure [5}. 

We searched for the remaining 64 spectra in the NED 
and foun d 14 of them in dif ferent publications, most o f 
them in IVigotti et all (|l997l ) and lLahulla et at! l|l99ll ). 
From these 14, only 4 could be identified as normal QSOs 
witho ut BALs: B3 0034+393 wit h z= 1.94 dOsmer et all 
1 19941 ); B3 0219+443 with z =0.85 jDiorgovski et al.lll990l ); 
B3 0249+383 with z= 112 dHenstock et al.lll997l ) and B3 
2351+456 with z=2.00 (|Stickel fc Kiihrlll993l ). The other 10 
are difficult to classify because they are either faint, with low 
signal-to-noise ratio in the continuum (just good enough for 
a correct redshift identification), or the wavelength cover- 
age does not include the region where the absorbing troughs 
could be present. Summing up, 59 out of 119 QSOs (roughly 
half of the sample) were inspected, with no one clearly show- 
ing the BAL phenom enon. These results j ustify the assump- 
tion that the list of l|Vigotti et al.|[l997l ) can be used as a 
representative sample of non-BAL quasars. 

Table [10] summarises the number of elements and the 
basic statistics on the spectral index for the samples un- 



der comparison. Si nce BAL QSOs are compact objects as 
iBecker et all (2000) discovered and we confirm in this work, 
we will also consider in our comparison a subsam ple of non- 
BAL QSOs extracted from IVigotti et~aLl (|l997t) which in- 
cludes only those compact sources with angular sizes <0.5 
arcsec at 1.4 GHz. The angular sizes at 1.4 GHz in Vigotti's 
sample were measured from maps taken w ith the VLA in C 
and D configurations l|Vigotti et al.|[l98^ 1 and for the most 
compact sources from observations in A configuration (priv. 
communication) . 

The first four rows in Table [TO] show that non-BAL 
QSOs have a median a?;| lower than BAL QSOs (—0.92 
and —0.53, respectively). The conclusion could be that BAL 
QSOs tend to have flatter spectra and are, on average, ori- 
ented at a lower angle with respect to the line of sight than 
normal QSOs. On the other hand, when the non-BAL QSOs 
sample is restricted to compact objects only, the median 
spectral index becomes —0.49, very similar to that of BAL 
QSOs. When the full sample of non-BAL QSOs is consid- 
ered, those extended objects have a higher contribution of 
emission produced in the lobes, which translates into steeper 
spectral indices for non-BAL QSOs. A similar behaviour can 
be found when analysing the rest-frame spectral indices afl 
and also ag 5 although in this last case the median values for 
BAL QSOs and compact non-BAL QSOs are more distant, 
i.e. —0.32 and —0.47 respectively. 

A Kolmogorov-Smirnov (K-S) test has been done in or- 
der to compare the af '\ spectral index distributions of these 
samples. The test yields a probability P<0.001 that both 
distributions are statistically equal, confirming that they are 
indeed different populations. A similar K-S t est re stricted 
to only compact sources from IVigotti et al.l l|l997l ) yields 
P=0.204 which is consistent with both samples having a 
similar distribution. The percentile plot in the upper panel 
of Figure [5] illustrates the shape of the spectral index distri- 
bution of these three samples. It can be seen that the im- 
portan t contribution of extended sources from lVigotti et all 
(1997) makes the percentile curve rise more rapidly than 
in the other samples. This result essentially means that the 
differences in spectral indices between BAL and non-BAL 
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Table 10. Sta tistics on spe ct ral in dices, BOO: iBecker et al.l 
fcOOCh: V97 : IVigotti et al.l l|l997t ): V97c: sources from 
IVigotti et al, | <1997ft with sizes <5 arcsec at 1.4 GHz 



Sctmplc 


N Min 


Max 


Mean 


Median 


Std 




Observer's frame: a [1.4 


- 8.4] GHz 




RBQ+BOO 


38 —1.50 


0.70 


—0.50 


—0.53 


0.43 


RBQ 


15 —1.19 


0.20 


—0.52 


—0.44 


0.34 


V97 


119 -1.47 


0.86 


-0.80 


-0.92 


0.43 


V97c 


50 -1.41 


0.86 


-0.53 


-0.49 


0.49 




Rest frame: 


a [6.0 - 


25] GHz 






RBQ 


15 -0.87 


0.17 


-0.38 


-0.32 


0.30 


V97 


51 -1.27 


0.17 


-0.71 


-0.82 


0.38 


V97c 


27 -1.21 


0.17 


-0.53 


-0.47 


0.40 




Rest frame: 


a [12 - 


25] GHz 
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15 -1.00 


0.45 


-0.40 


-0.58 


0.43 


V97 
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0.39 


-0.71 


-0.76 


0.42 


V97c 
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-0.62 
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Figure 9. Percentile plot of the samples involved in the K-S tests, 
illustrating the shape of the spectral index distributions. 



QSOs arise only because BAL QSOs are compact sources, 
while normal QSOs include both compact and extended. 

We will also compare the rest-frame spectral indices in- 
dices to look for possible intrinsic variations between BAL 
and non-BAL QSOs, once the redshift effect has been cor- 
rected for. The most reasonable comparison is using afl 
where the effect of synchrotron self-absorption seem less pro- 
nounced. Two K-S tests have been done comparing these 
samples. The first compares Vigotti's sources and RBQ 
yielding P=0.048 . while the second co mpares only the com- 
pact sources in IVigotti et all \ 19971 ) and RBQ, yielding 
P=0.629. Again we can talk about statistically different pop- 
ulations when all QSOs are included in the comparison, but 
not when only spectral indices of compact sources are com- 
pared. The percentile plot is shown in the lower panel of 
Figure [9] 

Even if the previous statistical test are consistent with 
similar spectral index distributions for BAL QSOs and non- 
BAL compact sources, it is interesting to note in Figure [9] 
that there might be a small deficit of steep sources among 
BAL QSOs with s pectral index ~ —0 .7, compared to the 
compact sources in IVigotti et all (Il997f l. This effect is more 
marked in the upper panel, but also insinuated in the lower 
panel where, nevertheless, the number of objects is signifi- 
cantly smaller. 

4.5 Polarisation 

The RBQ sample was searched in the NVSS database look- 
ing for polarised emission and only one BAL QSO, 1159+01, 
shows a high amount of linear polarisation at 1.4 GHz. For 
the rest of the sources only upper limits could be established 
typically below 1.5 per cent at this frequency. 

At higher frequencies and up to 43 GHz we have mea- 
sured the Sq and Su parameters as explained in Section 3, in 
order to obtain the degree of linear polarisation, m, and the 
polarisation angle X- Table [TT] shows the fractional polarised 
intensity, mostly upper limits, for the BAL QSOs in the 
RBQ sample. The only sources which are clearly polarised 
showing significant m at various frequencies are 1159+01 
and 1624+37. These sources are presented in different ta- 
bles (Table 1121 in this paper and Table 2 from iBenn et al.l 
2005). These two are the only bright sources displaying rel- 
atively high (i.e., >10 percent) polarised intensity at some 
frequency, having 1159+01 mi. 4=15 per cent and 1624+37 
with mio.5 = ll per cent. Apart from these, also the weak 
source 1625+48 with only 1.7 mjy at 43 GHz shows ni43 ~18 
per cent but this measurement is based on a > 3-<j detec- 
tion in Su and the uncertainty in this measurement is high. 
From the remaining sources in Table [TT1 1312+23 shows sig- 
nificant polarisation at 8.4 and 22 GHz, and there are 5 more 
sources with significant detections at only one frequency, in 
all cases below 4 per cent. 

It should be noted from Table[TT]that the most stringent 
upper limits are given by the VLA observations at 8.4 GHz. 
At this frequency the sources are bright enough to have reli- 
able detections and the relatively low noise achieved makes 
possible the determination of reasonable upper limits for the 
non-detections. At higher frequencies this is more difficult 
because the total power spectra drops down to fainter levels. 

At 8.4 GHz, one third of the sample (i.e, 5 out of 15) 
show at least a 3-cr detection in Sq, Su or both. From these 
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Table 11. Degree of polarisation, m„ (in percentage), at several frequencies, v (in GHz), for the RBQ sample. Most of the values are 3-<r 
upper limits and those showing a measured value can have a relatively large error but come from an at least 3-sigma det ection in Stokes Q 
or Stokes U. The BAL QSOs 1159+01 and 1624+37 presented separately (see Table [12] in this paper and Table 2 from lBenn et al.ll2005h 
are the only two sources with significantly detected polarisation. 



ID 




"12.6 


"14.8 


"18.3 


"18.4 


"110.5 


TO15 


"122 


"143 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


0039-00 


< 1.59 


~~ 


< 3.8 


< 14.5 


< 0.67 


<24.8 


<10.5 


<14.7 


< 8.3 


0135-02 


< 1.26 




< 3.1 


< 5.0 


< 0.87 


< 8.6 


3.9+4.0 


< 1.7 


< 7.4 


0256-01 


< 1.41 




< 4.5 


< 26.4 


< 3.9 


<16.4 


< 1.7 


<22.7 


< 7.2 


0728+40 


< 1.17 





<17.5 


< 30.7 


< 10.1 




<10.9 


< 5.1 




0837+36 




< 27.9 


< 2.2 


< 6.9 


< 0.6 


<16.4 


<27.8 


<15.2 


<10.1 


0957+23 


< 1.23 




< 0.7 


< 1.8 


< 0.5 


< 3.4 


< 0.9 


2.0+1.4 


< 5.6 


1053-00 


< 1.29 




< 4.2 


< 18.1 


1.3+1.0 


<20.8 


<12.8 


< 7.3 


<22.1 


1213+01 


< 1.29 




< 5.4 


< 13.7 


< 0.5 


<37.1 


<19.0 


< 9.9 


< 7.3 


1228-01 


< 1.35 




< 4.3 


< 10.7 


< 0.5 


<24.6 


<10.9 


< 1.0 


<23.7 


1312+23 


< 1.17 




< 2.5 


< 8.2 


2.7+0.7 


<16.2 


< 1.5 


3.6+2.3 


< 6.5 


1413+42 


< 1.26 




< 6.7 


< 8.9 


< 0.5 


<15.5 


<10.2 


< 1.1 


< 9.5 


1603+30 


< 1.32 


< 19.3 


< 2.4 


< 5.8 


1.0+1.0 


<14.0 


<11.2 


< 0.8 


<19.3 


1625+48 


< 1.08 


< 24.0 


< 6.4 


< 7.4 


< 4.9 


<16.6 


<19.9 


<19.8 


17.6+14.0 



a Degree of linear polarisation at 1.4 GHz extracted from the NRAO VLA Sky Survey, NVSS l|Condon et al.1 1998~) 



5 BAL QSOs only 1624+37 is strongly po larised showing 6. 5 
per cent of linearly polarised flux density dBenn et al.l l2005). 
The fractional polarised intensity in the other four objects 
(1053-00, 1159+01, 1312+23 and 1603+30) is below 3 per 
cent. The upper limits for most of the 10 undetected sources 
are below 1 per cent indicating strong depolarisation. 

For 1159+01 and 1624+37 it is possible the determi- 
nation of the Rotation Measure, RM, (x = Xo + RM ■ A 2 ) 
fitting a slope to the polarisation angles Xi as a func t ion o f 
the square of the observed wavelenghts. iBenn et al.l (|2005l ) 
already reported for 1624+37 a Rotation Measure of —990 
± 30 rad m -2 , which in the rest-frame translates into the 
extremely high intrinsic Rotation Measure, RM = -18350 
± 570 rad m~ 2 . For 1159+01 we have fitted the angles Xi 
in Table [T2j Our observations give a poor upper-limit at 
2.65 GHz, but an add itional measurement at th i s frequ ency 
has been obtained by ISimard-Normandin et al.l l|l98ll ). We 
get a best fit with a Rotation Measure RM=(-72.1 ± 1.4) 
rad m -2 and an intrinsic polarisation angle xo = —24 ± 
3 deg. This observed value should be corrected by the RM 
introduced by our own Galaxy, which is quite difficult to es- 
timate. To have a roug h estimate we have inspected the all- 
sky map presented by IWielebinski fc Krausd (|l993l ) where 
the RM of 976 extragalactic sources are plotted. The neigh- 
bouring sources around the position of 1159+01 have small 
(i.e., |RM| <30 rad m~ 2 ) Rotation Measures, which sug- 
gests that the Galactic contribution might be smaller than 
the determined value. Assuming no Galactic correction, the 
fitted Rotation Measure brought to the rest-frame would be 
higher by a factor (1+z) 2 becoming RM=(644 ± 12) rad 
m" 2 . 



5 DISCUSSION 

The morphologies and dimensions of the radio sources in 
BAL QSOs seem to be the keys to understand the orienta- 
tion and evolutionary status of these sources. There is ev- 
idence that BAL QSOs are associated with compact radio 
sources, which are supposed to constitute a high fraction 



Table 12. Polarisation properties of BAL QSO 
1159+ 01. The second measurement a t 2.65 GHz comes 
from lSimard-Normandin et al. 1 1119811 ) 



Frequency 


Telescope 


m 


X 


(GHz) 




(per cent) 


(deg) 


1.40 


NVSS 


15.03 ± 0.44 


-31.1 ± 0.6 


2.65 


Effelsberg 


<16 




2.65 




3.1 ± 1.9 


111 ± 17 


4.85 


Effelsberg 


2.1 ± 0.2 


-42 ± 6 


8.35 


Effelsberg 


0.95 ± 0.05 


-29 ± 5 


8.45 


VLA 


0.9 ± 0.4 


-29 ± 7 


10.5 


Effelsberg 


<1.2 




15.0 


VLA 


<1.2 




22.5 


VLA 


<0.4 




43.5 


VLA 


<1.2 







RM = 


-72.1 ±1.4 


rad m 2 




Xo = 


-24 ±3 


deg 











3 8.01 0.00 0.01 0.02 0.03 0.04 0.05 

A 2 (m 2 ) 



Figure 10. Determination of the RM for 1159+01. 
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of the young population o f radio sources, like e.g., CSS or 
GPS sources. iBecker et af] (2000) noted that 90 per cent of 
their sample of radio-loud BAL QSOs (extracted from the 
FBQS survey) present point-like structure at the resolution 
of 5" in contrast to the diversity of both point-like and ex- 
tended sources they found within the whole population of 
FBQS quasars. In fact, only a few BAL QSOs are known t o 
have an extended FR II radio structure (|Gregg et al.ll200rj ). 
Our 22-GHz observations of the RBQ sample confirm this 
tendency, because most sources appear unresolved at all fre- 
quencies constraining their apparent dimensions up to <0.1 
arcsec. This translates into projected linear sizes (LS) < 1 
kpc which are typical of GPS/CSS sources. It could be ar- 
gued that an extended component with a steep spectrum 
might be present, observable only at lower frequencies, and 
thus being the source sizes larger than those found at high 
frequencies. This hypothesis is a priori no supported by the 
fact that most sources are still point-like at 1.4 GHz as can 
be seen in the FIRST maps, with the exception of 1053—00. 

The radio spectra in Figure [5] also display the typi- 
cal convex shape of GPS-like sources, in most cases with 
a peak at ~ 0.5—10 GHz suggested by at least one or two 
points, or by upper limits at MHz frequencies. One con- 
cern about this result is that FIRST or WENSS flux densi- 
ties were observed at different epochs and variability might 
be an issue, although noticeable flux density variations are 
likely to be stronger in the optically thin region of the spec- 
trum. However, relatively deep integrations at MHz frequen- 
cies would be desirable to better constrain the shape of the 
spectra at low frequency. In addition, a project with VLBI 
mult i- frequency observations and the analysis on the pc- 
scale properties of several BAL QSOs of the RBQ sample 
has been started by us, and results will be presented in a 
subsequent paper. 

If radio sources associated to BAL QSOs are actually 
CSS / GPS sources, they should follow the anticorrelation be- 
tween int rinsic turnover fre quency and projected linear size 
found bv lFanti et all (Il990l) for CSS sources. This was, for 
instance, quantified by O'Dea fc Ba um ( 1997) in a combined 
sampl e of CSS and GPS sour ces compiled by iFanti et alj 
l|l99d) and IStangheflrm] i|l992l ). respectively. This relation- 
ship was found to be valid for both galaxies and quasars. 
The anticorrelation suggests that the mechanism producing 
the turnover simply depends on the source size. Figure 1111 
shows the relationship between intrinsic v pea k and projected 
LS for the samples of CSS and CPS sources mentioned be- 
fore, and we have also plotted for comparison our sample 
of radio-loud BAL QSOs. The turnover frequencies for the 
BAL QSOs appear in Table [9] while projected linear sizes 
have been extracted from the deconvolved sizes in Table [7] 

A group of 5 BAL QSOs match quite well the distri- 
bution of CSS/GPS sources, and they all show simple con- 
vex spectra, i.e., 0728+40, 1213+01, 1312+23, 1624+37 and 
1625+48. The three additional BAL QSOs displaying sim- 
ple convex spectrum, 0039-00, 0256-01 and 0957+23 are 
unresolved in our observations and their position can only 
be constrained by upper limits. Nevertheless these 8 sources 
with simple complex spectrum are consistent with the cor- 
relation. 

From the remaining 7 objects, 5 are resolved in our 
observations (0135-02, 0837+36, 1228-01, 1413+42 and 
1603+30) and they are located above the cloud of CSS/GPS 



sources. All of them have in common a complex spectrum 
with indications of a double component. In these sources 
the turnover has been determined using the high frequency 
component but it is possible in the spectra of 1413+42 and 
1603+30 to fit a second peak at lower frequencies (at 0.7 and 
0.8 GHz respectively). A second peak at low frequency could 
also be fitted to 1159+01 but we will not do that because 
being this source unresolved we cannot give the exact lo- 
cation in the diagram. When considering the low frequency 
peaks of 1413+42 and 1603+30 in the diagram, these two 
sources also agree very well with the anticorrelation as shown 
in Figure 1111 This better agreement could mean that the 
high-frequency part of these spectra could be dominated by 
the emission coming from a compact and active (possibly 
beamed) region smaller than the entire source, like a hot- 
spot. Unfortunately we do not have enough data to check 
whether this could be the case also for 0135-02, 0837+36 
and 1228-01. Low-frequency data will be important to sam- 
ple this part of the spectra to search for a possible second 
peak. Again, VLBI observations are crucial to further inves- 
tigate this interpretation. 

In fact, it has been shown that some GPS sources as- 
sociated to optical QSOs can be just normal flat-spectrum 
quasars for which a jet/knot/hot-spot dominates the spec- 
trum, which then adopts the charac teristic convex shap e 
of a GP S source, as discussed e.g . by ISnellen" et all ll 19991) . 
Recently IStanghellini et all (|2005h and lOrienti et all (|2006r i 
have studied VLBI samples of GPS sources and High Fre- 
quency Peakers (HFP, 5 GHz < < 10 GHz) confirming 
this scenario. In fact, these "contaminant" quasars seem to 
be present, even in high percentages, in many well-defined 
samples of GPS sources (|Torniainen et ai]|2005h . Two char- 
acteristics of "genuine" young G PS rad i o sou rces are low 
polarisation and low variability (o 'Dealll998l) . while flat- 
spectrum quasars show higher polarisation and variability. 

iTorniainen et all |2005), who have used an extensive 
database of sources with multifrequency data from several 
long-term variability programs, define as "bona fide" GPS 
sources those with a maximum variability Var^s <3. From 
our comparison in Section U which is, of course, based on 
only two epochs we see that many sources do not show sig- 
nificant variations at 1.4 GHz, and 3 out of 5 inspected at 
8.4 GHz show significant but not very strong variations. We 
cannot thus exclude these sou rces from the candidate s to 
young sources. As also shown in ITorniainen et alj (120051 ) the 
maximum variations are found at high frequencies, so multi- 
frequency multi-epoch observations of the RBQ sample will 
be very interesting for a more precise variability study. 

Compact GPS sources are supposed to be completely 
embedded in their narrow line region. On this basis, 
(|Cotton et al.l [2003) proposed the existence of a typical 
frequency-dependent scale, below which the GPS source be- 
comes completely depolarise d, the so-called "Cott on effect" . 
IFanti et all (|2004T ) and later iRossetti et all (120081 ) have con- 
firmed this behaviour anal ysing the polarisa tion properties 
of the B3-VLA CSS sample (|Fanti et al.ll200ll ). and they pro- 
pose models to describe the depolarising screen. The BAL 
QSOs in the RBQ sample have probably physical sizes of < 
1 kpc and if Faraday effects are important they should be 
completely depolarised even up to relatively high frequen- 
cies. This is the case for all sources except 1159+01 at 1.4 
GHz, and also for about 70 per cent of the sample at 8.4 
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Figur e 11. Rest-fram e turnover frequency versus projected linear size relation for samples of CSS and GPS sources (adapted from 
lO'Dea fc Baumlll997t ). Galaxies are represented by the symbol 'x', quasars by the simbol '+' and the RBQ sample is plotted as 
solid circles. Upper limi ts for those BAL QSO s unresolved at 22 or 43 GHz are marked by triangles. The dashed line represent the 
anticorrelation found bv lO'Dea &; Bauml l|l997f ): log v pea k = —0.21—0.65 log LS. In those two resolved sources with two peaks in the 
spectrum (1413+43 and 1603+30) the grey points indicate the new location in the plot when the lower frequency peak is considered. 



GHz. At these two frequencies the upper limits obtained for 
the degree of linear polarisation are reasonably small. From 
Table ITTI only 1312+23 seems to be weakly polarised at 8.4 
and 22 GHz, but no polarisation was found at 15 GHz. 

The two exceptions are 1159+01 and 1624+37. The first 
source is weakly polarised at 8.4 GHz but it is worth to 
note that the degree of linear polarisation increases as fre- 
quency decreases, becoming strongly polarised at 1.4 GHz. 
The opposite is expected when a depolarising plasma curtain 
is present. The reason for this anomalous behaviour might 
be that the emission at low frequencies arises from a differ- 
ent component, as suggested by the complex radio spectrum 
(see Figure [5} . An hypothetical extended structure with a 
relatively ordered magnetic field would contribute with sub- 
stantial linear polarisation, and this may only be detected 
at low frequencies if a steep spectral index is assumed for 
the emission of this extended component. 

The moderate observed Rotation Measure of (—72.1 
± 1.4) rad m -2 in 1159+01 is consistent with the fact 
that it is strongly polarised in the NVSS. Should the Ro- 
tation Measure be much higher, the large bandwidth used 
for the NVSS observations would produce a strong band- 
width depolarisation effect, probably hiding any measur- 
able polarised intensity. The extre me Rotation Mea sure of 
J1624+37, RM=-18350 rad m -2 (|Benn et al.ll2005l ) which 



is the second highest value known among QSOs and the 
highest among BAL QSOs, seems to be an exception and 
not representative of the BAL QSOs family. 

The median 777,8.4 of the 5 polarised BAL QSOs in the 
RBQ sample is ~1.3 per cent, but the remaining 1 BAL 
QSOs show m8.4 <f. As a comparison ISaikia et all (|l987f) 
found a median polarisation degree at 6 GHz of about 2 
per cent in a subsample of quas ars extracted fr om the sam- 
ple of ~400 compact sources of lPerlevI l|l982l ). ISaikia et all 
( 1987) found no significant differences between the median 
value of fiat-spectrum cores and CSS quasars whereas that 
radio sources associated to galaxies or to empty fields (no 
optical identifications) were found to be less polarised at 6 
GHz, with a media n value of ~0.5 per cent. More recently, 
IStanghellinil (|2003l ) found a mean fractional polarisation of 
1.2 and 1.8 per cent for GPS and flat-spectrum quasars, 
respectively, and m < 0.3 per cent for galaxies. 

The sources showing some polarisation in the RBQ sam- 
ple have probably polarised intensities consistent with flat- 
spectrum quasars. It is also likely that some of them (e.g., 
1159+01 or 1603+30) can show polarised emission due to 
Doppler boosted knots in jets, which is again compatible 
with their multi-component radio spectra peaking at high 
frequencies. However, our whole sample seem to be in bet- 
ter agreement with the GPS class than with the class of flat- 
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spectrum quasars, because of the high percentage of unpo- 
larised or weakly polarised sources. Of course, to firmly test 
this hypothesis it will be particularly interesting to do po- 
larimetric multi-frequency observations of a radio brighter 
sample, which is now available after the recent releases of 
SDSS. 

As far as the spectral index distribution is concerned, 
the different statistical tests suggest that the spectral index 
distributions of BAL and non-BAL QSOs are different. This 
difference is however due to the fact that all BAL QSOs in 
the sample have compact structures, while non-BAL QSOs 
have both extended and compact morphologies. When com- 
paring the samples of BAL QSOs and compact non-BAL 
QSOs, no significant differences are found in the spectral in- 
dex distributions. Thus, from our comparison no preferred 
orientation can be attributed to BAL QSOs. When larger 
samples become available, the comparison of spectral indices 
of BAL and non-BAL QSOs as a function of radio power 
and UV luminosity will be possible, allowing us to look for 
particular orientati ons in differ ent bins of luminosity, as has 
been suggested bv lElvisI (|2000f ) . 



6 CONCLUSIONS 

• A sample of 15 BAL QSOs associated with FIRST 
sources brighter than 15 mjy has been built. Radio contin- 
uum observations for these sources have been collected from 
2.6 up to 43 GHz in full polarisation. Flux densities cover- 
ing this radio range have been presented complemented by 
archive data at lower frequencies from 74 MHz up to 1.4 
GHz. 

• VLA maps in the most extended configuration show 
very compact morphologies for most sources at all frequen- 
cies being unresolved or slightly resolved at 22 GHz with 
a resolution of 80 mas. These translate into projected lin- 
ear sizes < 1 kpc, which are the typical sizes of CSS/GPS 
sources. 

• The spectra of these sources are typically convex, i.e., 
they seem to become flat or inverted at MHz frequencies 
probably due to synchrotron self-absorption, while at fre- 
quencies higher than 15 GHz they are definitely steeper. 
The spectra typically peak between 1 and 5 GHz in the ob- 
server's frame and about 1/3 of the sample shows complex 
spectra suggesting different components. 

• All BAL QSOs in our sample presenting simple con- 
vex spectra follow the anticorrelation between projected lin- 
ear size and turnover frequency found for CSS and GPS 
sources. Two BAL QSOs with complex radio spectra are 
also in agreement with this relationship if we consider the 
lower frequency peak of their spectra. The higher frequency 
peak might be in this two sources due to emission from a 
hot-spot or a knot in a jet. 

• A high percentage of the sample does not show sig- 
nificant variability when comparing the flux densities in 2 
different epochs. Some BAL QSOs show significant variabil- 
ity at 8.4 GHz but not strong enough to exclude them from 
the candidates to young radio sources. 

• Most BAL QSOs in our list are not strongly polarized 
either at 1.4 GHz or at 8.4 GHz, and sensible upper limits are 
given at these frequencies. Only two sources show significant 
polarisation at several frequencies. From these two, only the 



unusual BAL 1624+37 has a high Rotation Measures sug- 
gesting strong depolarization. The median fractional polar- 
isation of the sample is in better agreement to mean val- 
ues found in CSS/GPS sources, than those found for flat- 
spectrum quasars. 

• A series of statistical tests have been done to com- 
pare the spectral index distribution of BAL and non-BAL 
QSOs finding that these distributions are different. How- 
ever, no significant differences were found when comparing 
to compact non-BAL QSOs only. This is consistent with 
BAL QSOs spanning the same range of orientations as nor- 
mal quasars with respect to our line of sight to them. 
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